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ABSTRACT 

The masses of supermassive black holes (SBHs) show correlations with bulge properties in 
> ■ disk and elliptical galaxies. We study the formation of galactic structure within flat-core triaxial 

haloes and show that these correlations can be understood within the framework of a baryonic 
component modifying the orbital structure in the underlying potential. In particular, we find that 
terminal properties of bulges and their central SBHs are constrained by the destruction of box 
orbits in the harmonic cores of dark haloes and the emergence of progressively less eccentric loop 
orbits there. SBH masses, M., should exhibit a tighter correlation with bulge velocity dispersions, 
o"b, than with bulge masses, Mb, in accord with observations, if there is a significant scatter in the 
Mjj — oh relation for the halo. In the context of this model the observed M, — ctb relation implies 
that haloes should exhibit a Faber-Jackson type relationship between their masses and velocity 
dispersions. The most important prediction of our model is that halo properties determine the 
bulge and SBH parameters. The model also has important implications for galactic morphology 
and the process of disk formation. 

Subject headings: galaxies: evolution - galaxies: ISM - galaxies: kinematics & dynamics - galaxies: 
structure - hydrodynamics 

1. Introduction 

J JILA Visiting Fellow. Permanent address: Department 
of Physics & Astronomy, University of Kentucky, Lexing- The possibility that Supermassive black holes 



ton, KY 40506-0055, USA (SBHs) inhabit the centers of many if not most 

galaxies, and the observed correlation between 
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SBH masses and galactic bulge properties, 2 has 
potentially a fundamental significance for our un- 
derstanding of galaxy formation and evolution. 
The relationships between black hole and bulge 
properties include a loose relationship between 
SBH and bulge masses, M, <~ 0.001Mb, and an 
apparently much tighter one between the SBH 
mass and the velocity dispersion in the corre- 
sponding bulge, M, <~ cr B (e.g., Ferrarese & Mer- 
ritt 2000; Gebhardt et al. 2000; Tremaine et al. 
2002; cf. reviews by Kormendy & Gebhardt 2001; 
Merritt & Ferrarese 2001). 

In this paper we attempt to provide a physical 
explanation for these relationships between SBHs 
and their host galaxies. Our model is based on 
the interaction between the dark haloes of galax- 
ies and the baryonic components settling in their 
midst. As baryonic matter accumulates to form 
the bulge and SBH, the orbital structure of the un- 
derlying gravitational potential is modified. This, 
in turn, affects the subsequent accumulation of 
gas, which is highly dissipative and therefore is 
sensitive to orbital geometry. Consequently, halo 
properties will determine those of the bulge and 
the SBH, and this will give rise to correlations be- 
tween them that are compatible with the observed 
ones. 

Cuspy, triaxial haloes appear to be a natural 
outcome of dissipationless collapse in a cold dark 
matter (CDM) dominated universe (e.g., Cole & 
Lacey 1996). Interactions with the baryonic com- 
ponent during the initial stages of collapse may 
affect the triaxiality, making it milder but still 
non-negligible (Dubinsky 1994). Before it even be- 
comes dominant in the central regions, a clumpy 
baryonic component can also level off the cen- 
tral cusps of dark haloes, producing harmonic 
cores (El-Zant, Shlosman & Hoffman 2001). Even 
though this may affect the equidensity contours, 
again making them rounder, it need not sym- 
metrize the equipotentials; these can remain asym- 
metric if triaxiality is not affected beyond some ra- 
dius (for example, a homogeneous bar has a non- 
axisymmetric force contribution inside its density 
figure). From an observational standpoint it ap- 
pears that haloes of fully formed galaxies tend 
to have nearly constant density cores (de Blok & 

2 The term "bulge" in this paper will refer to bulges of disk 
galaxies, or to elliptical galaxies, when no disk is present. 



Bosma 2002) and that residual potential axial ra- 
tios of about 0.9 in CDM haloes are plausible, even 
in -present day galaxies (e.g., Kuijken & Tremaine 
1994; Rix & Zaritsky 1995; Rix 1995). 

The orbital structure of the inner regions of 
slowly rotating, non-axisymmetric potentials with 
harmonic cores is dominated by box orbits (e.g., 
Binney & Tremaine 1987), which have no particu- 
lar sense of circulation. They are self-intersecting 
and, therefore, cannot be populated with gas. Dis- 
sipation causes material to sink quickly toward the 
only long-lived attractor available — the center 
(e.g., Pfenniger & Norman 1990; El-Zant 1999). 
In the process, interaction with the triaxial har- 
monic core causes the baryonic material to lose 
most of its angular momentum. The combined 
system would, therefore, also be slowly rotating. 
Thus, unless star formation terminates the col- 
lapse, the final concentration of the first baryonic 
material could be extremely large. 

The onset of star formation is expected to 
occur when the baryonic material becomes self- 
gravitating — roughly speaking, when its density 
becomes larger than that of the halo core. This 
also happens to be the criterion for the destruction 
of the harmonic core and the emergence of loop (or 
tube) orbits, which do have a definite sense of cir- 
culation. The role of the SBH is to contribute to 
the emergence of these orbits in the very central 
region. It is the collusion between the SBH and 
the more extended hot baryonic component in cre- 
ating the loop orbits that leads to the correlations 
claimed in this paper. 

Whereas, in three dimensions, a box orbit can 
be represented superposition of three (gen- 
erally) incommensurable radial oscillations along 
mutually perpendicular axes, which are efficiently 
attenuated by dissipation, loop orbits are best de- 
scribed in terms of modest radial and vertical ex- 
cursions superposed on rotational motion about 
the center. While the vertical and radial ex- 
cursions, like the oscillations characterizing box 
orbits, are attenuated by dissipation, the rota- 
tional motion is not efficiently dissipated among 
gas clouds populating such orbits (in the same 
sense of circulation) . 

The minimization of the radial and vertical os- 
cillations results in closed periodic orbits that are 
confined to a plane. We expect the amount of gas 
dissipation on loop orbits to depend on their ax- 
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ial ratios — for, again, motion on highly eccentric 
orbits would lead to shocks and the accompanying 
loss of angular momentum on a short dynamical 
timescale, as observed in numerical simulations of 
gas flows in barred galaxies (e.g., Heller & Shlos- 
man 1994). Lacking a detailed model for the dis- 
sipation rate associated with gaseous motion, we 
will assume that there is a critical eccentricity be- 
low which the loop orbits can serve as long-lived 
attractors for dissipative motion. In other words, 
gas populating these orbits will evolve only secu- 
larly and not dynamically. Fortunately, this de- 
pendence on critical eccentricity will turn out to 
be weak. 

In the following we explore, within the above 
framework, the formation of galactic bulges and 
central compact objects. We show that it is possi- 
ble to deduce a well-defined linear correlation be- 
tween the masses of the central SBHs and corre- 
sponding bulges in given dark matter haloes, with 
mass ratios comparable to those observed. We 
demonstrate that, if dark matter halo cores follow 
a Faber- Jackson type relation (Faber & Jackson 
1976) between their masses and velocity disper- 
sions, then a similar relation also applies to the 
bulges. An M, — ctb relation between the SBH 
mass and the bulge velocity dispersion, which un- 
der certain assumptions can have smaller scatter 
than the M, — Mb relation, also follows naturally. 
Within this framework, it is possible to make a 
number of testable predictions concerning the re- 
lated structures of bulges, SBHs and their host 
haloes. We discuss these in the final section. For- 
mal aspects of the perturbation analysis applied 
to bulgc-halo systems have been deferred to Ap- 
pendices A and B. Preliminary results of this work 
have been reported by Shlosman (2002). 

2. The model 

Since, at this stage, we are interested in generic 
dynamical phenomena related to orbital structure, 
the exact form of the halo potential is immaterial 
- as long as it exhibits a harmonic core where 
no loop orbits can exist. The exact distribution of 
the baryonic material is also not crucial, except, 
again, for its central density distribution. If this 
diverges, say as p oc r -1 , loop orbits will be cre- 
ated all the way to the center, since the bulgc-halo 
system no longer possesses a harmonic core. In 



this case there is no need to form the SBH. If, on 
the other hand, the (proto)-bulge has a harmonic 
core of its own, there will still be a nearly constant 
density region near the center with only box orbits 
— unless a central point mass is present. We will 
assume that 'cuspy' bulges, when they exist, are 
products of later evolution (for example, a result 
of star formation and cold dissipationless collapse: 
e.g., Lokas & Hoffman 2000). 

We use a logarithmic standard form (e.g., Bin- 
ney & Tremaine 1987) to represent the halo po- 
tential: 

$ H = \vi io g (i?| + .x 2 + /r V + 7- V), (i) 

where Vh is the asymptotic (in the limit R 3> i?H 
and /3, 7 — > 1) circular velocity, i?n is the core ra- 
dius and /3, 7 < 1 are the potential axis ratios. We 
will consider the process of bulge formation to be 
terminated, or at least substantially slowed down, 
when the combined (baryonic plus halo) potential 
admits sufficiently round non-intersecting loop or- 
bits, which permit the long-term circulation of gas 
without excessive dissipation. Such motions neces- 
sarily take place in a symmetry plane determined 
by the (orbit-averaged) angular momentum (e.g., 
Frank, King & Raine 2002). Therefore, for our 
purposes, it will suffice to consider only orbits in 
the plane z = and to ignore the vertical dimen- 
sion. We fix (3 at 0.9 (though the effects of its 
variation are discussed where relevant) and adopt 
a threshold axial ratio p cr ;t for orbits that can be 
populated with gas. The value of p cr j t that best 
describes when this happens remains to be investi- 
gated, but its mere existence is what is important 
here. As we show below, our results are rather 
insensitive to the exact value of p C rit- 

The halo core mass will be taken to be Mh = 
V^Rh/G and we define the halo density to be 
Ph = M H /i?H = V^/R^G. This is conveniently 
close to the value of the density in the region 
where the potential is effectively harmonic — that 
is, within the region where, in the absence of the 
bulge component, no loop orbits exist for (3 = 0.9. 
We note, however, that a bulge with a larger core 
radius probes a region of the halo core with smaller 
average density than the region probed by a bulge 
with a relatively small core radius. This effect will 
be discussed in Section 4. 

Since our aforementioned criterion depends on 
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the potential in a chosen symmetry plane, the 
exact three-dimensional mass distribution of the 
bulge is unimportant, as the same planar potential 
may arise from a variety of these. A particularly 
convenient form for the potential in this situation 
is that of Miyamoto & Nagai (1975), 



\jx 2 + y 2 + (A + VB 2 + z 2 ) 

This potential can approximate a disk-bulge sys- 
tem, with the parameters A and B determin- 
ing the scale length and height, respectively. In 
the symmetry plane, therefore, the potential can 
arise from a range of density distributions — from 
highly flattened to spherical ones. In general, 
when we refer to the "bulge" we will have a spher- 
ical system in mind with A = 0. In particular, all 
densities quoted are calculated under this assump- 
tion. We will comment on the effects of assuming 
a flattened distribution where relevant. In that 
case Md in eq. (2) is replaced by Mb, the spher- 
ical bulge mass. The bulge core is defined by the 
radius Rb = A + B and encloses the core mass 
M BC = 2- 3 / 2 M B . The bulge core density is de- 
fined as pB = Mbc/R b - 

The scaling relations that arise do not depend 
on the absolute masses and lengths, but instead, 
on the ratios of these quantities among the vari- 
ous contributions to the potential. An important 
property of the dynamics in our model is that any 
transformation that leaves the relative masses and 
lengths of the components constant will also not 
affect the orbital structure once scale transforma- 
tions are taken into account. Such transforma- 
tions only change the time units, hence the dy- 
namics remains invariant on a new characteristic 
timcscale. In addition, in regions where the bary- 
onic and halo components have the same mass dis- 
tributions, i.e., inside their cores, only the density 
ratios of these components determine their rela- 
tive contributions to the potential, and hence the 
orbital structure. 

For the remainder of this paper, we will de- 
scribe all components of the potential in terms of 
units scaled to the dimensions of the halo poten- 
tial, cq. (1). We scale all distances to i?n, with 
R B = Rb/Rh, and all masses to M H so that 
M B = Mb I Mb- The normalized bulge core den- 
sity is then p B = 2~ 3 / 2 M B /-R B . The bulge poten- 



tial in scaled units is 

*b = -Vj / 23/2fa ^ B ■ (3) 
y/x* + y 2 + r 2 b 

The black hole potential, in the same units, is 

.. = -vi^sA, <4) 

v x + y 

where K = M./M B - Note that all orbits in the 
total potential <J> = $h + + depend only on 
three parameters, which may be taken to be Rb, 
pB, and K. 

3. Bulge formation in triaxial haloes and 
the M, — Mb correlation 

In this section, we demonstrate that the argu- 
ments outlined above imply a linear relation be- 
tween the SBH and bulge masses, provided that 
the baryonic and CDM density distributions are 
not completely axisymmetric and exhibit, at some 
stage, nearly constant density cores. 

3.1. Critical eccentricities and the emer- 
gence of SBH-bulge correlations 

3.1.1. The inner minimum and the role of bulge 
density 

Fig. 1 (the result of orbital integrations) ex- 
hibits the axis ratios of the closed loop orbits 
P = £max/j/max as a function of the longer axis 
length (t/max)- In this figure, we have fixed the 
scaled bulge core density, pB = 5.52, and the SBH- 
to-bulgc mass ratio, K = M./M B = 10~ 3 , but 
vary the scaled bulge core radius, Rb- Conse- 
quently, the bulge masses Mb range over a fac- 
tor of 100. As expected, when only the halo con- 
tributes to the potential (dotted line), p quickly 
drops to zero inside the halo core, since no loop or- 
bits can exist there. (For halo potential axis ratios 
smaller than the adopted value of 0.9, i.e., greater 
triaxiality, p reaches at radii that are progres- 
sively closer to the halo core radius, y m ax — > 1. As 
the baryonic component, both in the bulge and the 
SBH, becomes progressively more massive, corre- 
spondingly rounder loop orbits appear. Loop or- 
bits in the very central region are produced by the 
presence of the SBH. The latter is represented nu- 
merically by a Plummer sphere with a softening 
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scale of 0.002, in scaled units. Near the center, 
where the SBH contribution is dominant, the axis 
ratios of the loop orbits tend to p — > 1. 
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Fig. 1. — Axis ratios of closed loop orbits p = 
imax/j/mra as a function of ymax- The dotted line 
corresponds to a system with only a halo contri- 
bution (i.e., no bulge and no SBH), the solid line 
refers to a system with a "bulge" in the form of 
a Plummer sphere with Mb = 1 and Rb = 0.4, 
the dashed line represents a bulge with Mb = 0.1 
and i?B = 0.19, while the dashed-dotted line rep- 
resents a bulge with Mb = 0.01 and a radius of 
.Rb = 0.086. In all cases with the bulge present, 
Pb = 5.52 and an SBH with a fraction K = 10~ 3 
of the bulge mass is present. 

The axis ratio curves exhibit two distinct min- 
ima. It is apparent that the inner minima of the 
curves in Fig. 1 have (nearly) the same values of 
p, but occur at different radii oc Rb- This is a 
consequence of assuming a fixed ratio of SBH-to- 
bulge mass, K, and a constant bulge core density 
Pb > I- Because the bulge and halo densities are 
nearly uniform in the inner regions, their relative 
contributions to the potential are determined by 
their density ratio. 

In the absence of the SBH contribution, p in 
Fig. 1 would tend to zero within the effective har- 
monic core of the bulge- halo system (see cq. [B8]). 
When the relative density is kept constant (and 

~ 1/3 

> 1) this is proportional to Rb ~ M B ' (the ef- 
fect of density variation will be examined in Sec- 
tion 3.2). The reason for this is that inside the 



nearly constant-density bulge core, the potential, 
which is now a superposition of two nearly har- 
monic potentials, is nearly harmonic (even if the 
new core is less triaxial than the halo alone, e.g., 
if the bulge is assumed to be spherical) . The inner 
minima of the curves in Fig. 1, therefore, corre- 
spond to a transition from the region where the 
SBH provides the dominant contribution toward 
the creation of loop orbits to that where the bulge 
provides this contribution. Therefore, the min- 
ima occur at radii where the gravitational acceler- 
ation due to the SBH is proportional to that due to 
the bulge+halo, i?. - (M./M B ) 1/3 i? B = i^ 1/3 i? B 
(since the minimum is located well within the 
bulge core radius). Since K is taken to be con- 
stant, we expect the minima to occur at radii 
oc Rb- 

In Section 2, we have defined the critical value 
Pcrit above which gas circulation can be sustained 
for secular timescales. Here we have shown that 
given a critical value, p C rit, for the inner mini- 
mum and the relative bulge-to-halo core density 
ratio, a value of K associated with this minimum 
is determined. However, satisfying the condition 
P > Pcrit at the inner minimum does not imply 
that the same condition is satisfied at all radii 
within Rb- Until it is satisfied at the outer mini- 
mum, as well, further evolution can occur. This is 
briefly sketched below and outlined in more detail 
in Section 3.4. 

3.1.2. Outer minimum, minimal bulge mass, and 
final black hole mass 

Until the critical value of p is reached at all 
radii inside the halo core, the bulge will continue 
to grow, since outside the bulge core p declines 
again. It is clear from Fig. 1 that, in contrast to 
that of the inner minimum, the value of p at the 
outer minimum depends more sensitively on the 
bulge mass than on its density. The growth of 
the bulge, therefore, determines a minimal bulge 
mass fixed by the condition that p > p cr ; t at all 
radii within the halo core. 

If the SBH did not continue to grow along with 
the bulge, this growth in the bulge mass would 
decrease the linear correlation coefficient K, but, 
as we show in Section 3.4, this leads only to near- 
independence of the final value of K on jj cr ;t and 
still keeps the values of K within the observed 
range. 
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3.2. Scaling relationships determined by 
inner minima 

As shown in the previous section, the creation 
of loop orbits with axis ratio above a given value 
depends, in the central region, solely on the ratio 
of the bulge-to-halo core density and the mass ra- 
tio K = M./M B . The actual value of the bulge 
mass determines only the position of the minimum 
(as a fraction of the halo core radius), not the axis 
ratio at the minimum. But even the radius of the 
minimum is largely insensitive to the bulge mass, 

1 /3 

i.e., oc M B ' , provided that the minimum actually 
exists and the bulge core density relative to the 
halo harmonic core density remains constant. We 
will now examine the bounds on M, and K which 
ensure the existence of a minimum p > p cr it inside 
the bulge core and observe the effect of varying 
the bulge density. 

We calculate and plot K = K(pb) for different 
choices of p cr it (Fig. 2), in the following way. First, 
without any loss of generality, we choose a bulge 
mass Mb that is large enough that the outer min- 
imum satisfies p > p cr it, for all p cr ;t < 0.9. This is 
done in order to focus on the effects of the inner 
minima only. (Recall, from the previous section, 
that the value of the outer minimum depends on 
the mass of the bulge.) Fig. 2, therefore, exhibits 
the effect of pb and K on the inner minima of the p 
curves at R < Rb- In order to obtain K = K(pb), 
we vary M., which is some fraction K of the cho- 
sen bulge mass. The density is varied by contract- 
ing the bulge (by decreasing its core radius) until 
the condition p > p CI i t is satisfied everywhere in- 
side the bulge core. 

One observes in Fig. 2 that at smaller densi- 
ties K tends toward an asymptotic (and maximal) 
value associated with a given p CI1 t- We refer to 
this hereafter as the "asymptotic" regime. In the 
transition to this regime, along a p cr ;t = const, 
curve, the value of K is monotonically increas- 
ing. In other words, the bulge contribution to- 
ward the creation of loop orbits of a given eccen- 
tricity gradually diminishes, and is compensated 
by a greater contribution from the SBH compo- 
nent. In the process, the minimum in p moves 
outward. The asymptotic regime corresponds to 
a situation in which the existence of loop orbits 
with the required elongation depends on the value 
of K, irrespective of the bulge density (that is, 
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Fig. 2. — SBH-to- bulge mass ratio K for different 
Pcrit as a function of the bulge-to-halo core density 
ratio required to create closed loop orbits with axis 
ratios p > p cr it within the bulge core, R < Rb- Mb 
is chosen to satisfy p > p cr it everywhere outside 
the bulge. From left to right, the characteristic 
p cr it associated with a given curve increases from 
0.3 to 0.9, in increments of 0.1. The bulge core 
density pb is varied by fixing Mb and changing 
^b- The vertical dashed line at /3b = 1 is the 
approximate boundary between the self- and non- 
self-gravitating regimes in the bulge. For densities 
Pb ~ 1, the pcrit = 0.9 curve corresponds to an 
(unrealistically) overmassive SBH and is shown for 
comparison only. The range of K values, for these 
densities, thus is limited to within 10~ 4 — I0~ 2 . 

effectively, only on the SBH mass). In this limit 
the SBH contribution to the potential is sufficient 
to create loop orbits with p > p cr ;t at all radii 
within the halo core (i.e., at R < I), without ad- 
ditional contributions from the bulge component 
(cf. Fig. 3). Essentially, it corresponds to a rapid 
collapse to the center and formation of the SBH by 
a large fraction of the baryonic material, bypassing 
the formation of a bulge. This regime appears to 
be of academic interest only, as it implies pB < 1 
all the way to the center; it will not be discussed 
further. 

A second regime characterizes the dynamical 
state of the SBH-bulge-dark halo system with 
pB > 1. In this "scaling regime" the log if — 
log ,6b curves are parallel straight lines with av- 
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Fig. 3. — M. sufficient to create loop orbits 
rounder than a given p cr i t without collapse of the 
baryonic material to form a bulge. 

erage slopes of about —2.5. This can be ex- 
plained in the following manner. From eq. (B8) 
we know that the radius of the effective harmonic 
core of the bulge- halo system, y max , is propor- 
tional to fig^ifg 1 ' 2 . In the scaling regime, the 
location of the minimum in p will be proportional 
to y max . Now, in order to maintain the minimum 
at a specified value of p cr ; t , the gravitational ac- 
celeration due to the SBH must be proportional 
to the acceleration due to the halo (which con- 
tains the nonaxisymmetry) at y max . This implies 
PbRbK/v^ oc y max . Substituting for y max and 
noting that pb oc Mb/R b , we obtain K oc /5 B 5 ^ 2 . 
This is approximately what is found from Figure 2 
and it also holds if we change the mass of the 
bulge keeping the radius constant and again us- 
ing eq. (B8). Note, furthermore, that for constant 
Pb this relation predicts the radius of the effective 
harmonic core (and the location of the minimum) 
to be proportional to Rb, as expected from the 
heuristic considerations of the previous section. 

While the asymptotic values of K depend on 
Mb, the values of K in scaling regime do not. The 
transition between these two regimes can, there- 
fore, be characterized by a sharp change in the 
behavior of K, as seen from Fig. 2. 

We are mainly interested in the scaling re- 
gion, because the onset of star formation can be 
tied to the baryonic component becoming self- 



gravitating, which corresponds roughly to the 
bulge core density exceeding that of the back- 
ground halo, i.e., pa > 1. In this case, and for 
Pcrit 0.8, the SBH contributes significantly to 
the potential only at radii The p curves, 

therefore, exhibit a definite inner minimum well 
inside the halo core. 

The values of ps and the initial K, in princi- 
ple, also depend on the critical eccentricity, p C rit, 
of the inner minimum, which is expected to be 
independent of the bulge mass, but which does 
depend on complex gas dynamics. If, after the 
SBH forms, the bulge mass falls short of the value 
required for creating sufficently round closed loop 
orbits at all radii inside the halo core, it will con- 
tinue to grow until the outer minimum of the axis 
ratio curve also attains p > p C rit- In Section 3.4, 
we show that the SBH growth, if continued, be- 
comes intermittent. In the next section, we also 
demonstrate that one can place constraints on the 
possible range of K values by considering this con- 
straint on the outer minimum. 

3.3. Minimal bulge mass and bulge growth 

To obtain loop orbits rounder than a given p cri t 
at all radii R < 1, and not only in the central re- 
gion, one in fact needs to take into account the 
bulge-to-halo core mass ratio and not only the 
density ratio. This can already be seen in Fig. 1, 
where all curves with SBH and bulge contributions 
exhibit an inner minimum with p > 0.9, while p 
declines significantly as one moves further out. For 
low bulge masses, an additional outer minimum 
emerges, within Rb < R < 1, before p rises again 
outside the halo harmonic core. Only one of the p 
curves, corresponding to the most massive bulge, 
has p > 0.9 at all radii. Therefore, for a given 
halo, there is a minimal bulge mass, Memm, that 
is necessary to create sufficiently round loop orbits 
at all radii R > Rb- 

If we assume that the bulge core density varies 
at most by a factor of a few (/5b = 0[1]), Rb, 
which is only weakly dependent on the bulge mass 
and density (_Rb ~ M^ 3 Pg 1 ^ 3 ), varies little. The 
minimal bulge mass is also nearly constant (calcu- 
lations show, for example, that it varies by ~ 30% 
when pB changes from 1 to 2.5). There is more 
sensitivity to the assumed p cr it- In Fig. 4 (aster- 
isks) we show the bulge mass required to create 
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loop orbits with p > p cr it at all radii outside Rb , 
for pb = 1. These are the masses necessary, at 
this density, to produce outer minima with the re- 
quired values. For bulge densities pb ~ 1 these 
minima lie at radii R > 0.5. Thus, unless the 
SBH can contribute significantly to the potential 
at radii comparable to Rb, which seems implau- 
sible, the dominant contributions to the potential 
at the radii examined here should be only those 
due to the bulge and halo components. Therefore, 
once the halo parameters are fixed, the creation of 
loop orbits with given p > p cr it outside the bulge 
core will depend only on the bulge parameters Rb 
and Mb- 



contribute significantly at all radii — no matter 
how massive the bulge is. (For p B = 1, the curve 
with p cr it = 0.9 in Fig. 2 lies in this regime.) Fur- 
thermore, the minimal bulge density required to 
create loop orbits with p > p cr iti without appeal 
to an overmassive SBH, increases rapidly when 
Pcrit ^ 0.8 (Fig. 5), suggesting that the bulk of 
the material may form stars before the condition 
P > Pcrit is reached everywhere within the halo 
core. If we demand that pb < 1, then p cr it must 
be smaller than 0.8 in order for the model to be 
plausible, implying that K must lie in the range 
io- 4 - 10~ 2 . 



* * * M Bmir/ M h 
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Fig. 4. — Asterisks: bulge masses MBmin(pcrit) 
required to produce loop orbits with p > p cr it at 
radii outside the bulge core Rb, as a function of 
Pcrit for Pb = 1- Filled circles: corresponding ratio 
of the SBH-to-bulge mass K required to produce 
inner minima in the loop orbit axis ratios with p > 
Pcrit as a function of p crit at p B = 1 (from Fig. 2). 
The p cr it > 0.8 points are shown for comparison 
only, as this condition requires (for pb ~ 1) an 
unrealistically massive SBH, with strong effects on 
the potential at all radii within Rb- This accounts 
for the steep increase in K for these values of p cr it- 

For pb — 1, a well-defined inner minimum in 
the p curve exists inside the bulge core only if 
Pcrit < 0.8. For larger p cr it, the inner minimum 
moves outward and merges with the outer one. 
Systems with these properties lie outside the scal- 
ing region in Fig. 2. In this case, for the con- 
dition p > p cr it to be satisfied, the SBH has to 
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Fig. 5. — Minimum density required for creating 
loop orbits with p > p cr it without appeal to an 
overmassive SBH (see text). 

On the other hand, star formation can be some- 
what delayed to higher densities. A reasonable 
range lies within pB ~ 1 — 10. Still higher den- 
sities in the bulge can be excluded based on the 
observed rotation curves. In this range of pb , val- 
ues of Pcrit > 0.8 become feasible without invoking 
unrealistically massive SBHs, as seen from Fig. 5. 
One may also reasonably assume that, due to en- 
hanced dissipation, higher densities require larger 
Pcrit to maintain long-lived gaseous motion (de- 
tailed modeling, however, will be necessary to de- 
termine exactly how p cr it depends on the system 
parameters). In this situation, the range in K 
that can be inferred from Fig. 2, is again about 
10~ 4 — 10~ 2 and compatible with the observed 
value of <~ 10 ~ 3 , which exhibits a significant scat- 
ter. 
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The values of K obtained so far, which should 
be regarded as initial values, are fixed by the in- 
ner minimum of the axis ratio curves. As argued 
in Section 3.1, further evolution of the bulge will 
occur unless sufficiently round loop orbits also ex- 
ist at the outer minimum. In the next section we 
examine how this evolution can modify the range 
of K, if at all. 

3.4. Black hole growth and SBH-to-bulge 
mass ratio 

There are basic differences between the condi- 
tions for satisfying p > p cr i t at the inner and outer 
minima in Fig. 1. As discussed in Section 3.1, 
the initial growth of the SBH is terminated after 
the inner loop orbits reach the critical eccentricity. 
The growth of the bulge, on the other hand, con- 
tinues until the required p CI1 t is reached at all radii 
inside the halo core. Since the increase of bulge 
mass, at roughly constant density, moves the in- 
ner minimum outward as described in Section 3.1, 
the SBH would have to continue growing in pace 
with the bulge mass if it is to maintain the critical 
value of K (which is determined only by the values 
of pb andp cr it). Does this actually happen, and if 
so, how? 

To see how the parallel growth of the SBH 
and bulge could come about, consider the solid 
curve in Fig. 6, which represents the state of the 
SBH/bulge system at the end of the initial stage of 
infall. This is equivalent to one of the curves (say, 
the dash-dotted line) in Fig. 1, except that we have 
chosen p cr ;t ~ 0.6 for the inner minimum, instead 
of 0.9. This line has two maxima: a left-hand 
one at the origin, determined by the SBH, and 
a right-hand one, determined by the bulge core. 
Note that the outer part of the solid curve drops 
below p cr it, implying that the bulge will continue 
to grow. Now suppose this growth occurs without 
corresponding growth of the SBH (in contrast to 
the case in Fig. 1, where K is kept constant). If 
the bulge continues to grow while the SBH growth 
is stopped, the right-hand maximum moves fur- 
ther to the right, while the left-hand maximum 
stays the same. This opens a widening 'gap' at 
the position of the inner minimum which drops 
below p cr it, allowing gas to flow again toward the 
radius of influence of the SBH, R,. This is illus- 
trated by the dashed curve in Fig. 6, which cor- 
responds to a factor of 10 increase in Mb accom- 
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Fig. 6. — Axis ratios of closed loop orbits p = 
imax/i/max as a function of y max (as in Fig. 1 but 
for pb = 1 and growing the bulge mass without 
that of the SBH mass, thus decreasing K). The 
solid line corresponds to K = 10~ 3 and Mb = 
10~ 2 . The dashed line corresponds to K = 10~ 4 
and M B = 10 _1 . Note the 'gap' which opens with 
decreasing K and allows the gas to flow in again 
as p becomes less than p crit in this region. 

panied by a similar drop in K. The gas is ex- 
pected to accumulate in this vicinity. However, 
once a substantial amount of gas has collected 
at R,, it becomes self-gravitating and is prone to 
global self-gravitating instabilities. The fastest of 
these instabilities, m = 2 modes or bar instabil- 
ity (e.g., Bardeen 1975), have been discussed in 
a similar context by Shlosman, Frank & Begcl- 
man (1989). They induce rapid (dynamical) gas 
inflow. Hence we expect that the growth of the 
SBH at this stage will be intermittent, but because 
the time-averaged conditions for infall at the inner 
and outer maximum do not change, we expect the 
value of K to stay roughly constant within the 
range given in Fig. 2, namely ~ 10~ 4 — 10 -2 , de- 
pending on the value of p C rit- 

Is there always enough gas to accumulate at 
~ R t in order to trigger a bar instability, due to 
the opening of the gap between the SBH and the 
bulge? One can imagine the opposite extreme to 
that discussed above, in which insufficient gas en- 
ters the widening gap from outside, perhaps be- 
cause star formation is efficient during the early 
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stages of infall. Let us suppose, for the sake of 
argument, that the SBH does not grow beyond its 
initial mass, as determined by the initial value of 
K. How does this affect the range of final K val- 
ues? 

Figure 4 shows minimal bulge masses, A^Bmin ^ 
1CT 1 — 10 ~ 3 , for a range of p cr i t values. While 
we have no estimate for the initial bulge masses 
which define the initial K in Fig. 2, it is possible to 
rule out very small masses, e.g., < 10~ 3 , because 
the total baryonic mass within the halo core is 
taken to be about 10%. This means that for lower 
Pcrit ~ 0.3—0.4, the initial bulge mass is equal to or 
larger than Memm estimated in Section 3.3, and, 
therefore, will not grow beyond its initial value. 
It is easy to understand this result, because one 
needs small baryonic masses to create loop orbits 
with such large eccentricities. This means that 
in this regime the initial value of K ~ 10~ 4 (in 
Figures 2 and 4) is also its final value. 

Alternatively, in the regime of larger p CIlt w 
0.7 — 0.8, the bulge can grow at most by a factor of 
~ 30— 100, reducing the initial K by this amount. 
Luckily, the initial values of K for high p C rit lie 
around their high end, ~ 10 ~ 2 . Reduction by up 
to two orders of magnitude brings them again to 
about 10~ 4 . Hence, whether the SBH grows af- 
ter the initial stage or not does not destroy the 
M. — Mb correlation and does not move the val- 
ues of K outside the observed range. A corollary 
of this discussion is that the final value of K ap- 
pears insensitive to the value of p cr it, if the growth 
of the SBH is supressed, and depends strongly on 
p cr it, if the SBH grows in tandem with the bulge. 
Note also that the variation of MBmin(Pcrit) (rep- 
resented by the asterisks in Fig. 4), follows closely 
that of K (pcrit) (represented by filled circles). The 
ratio of these two quantities, M.MH/M| min <~ 
0.03 is, therefore, largely independent of p C rit- In 
other words, the minimal bulge mass is predicted 
to be about 5 times the geometric mean between 
the SBH and the halo masses, provided that the 
SBH grows along with the bulge. 

To summarize, the initial value of the black 
hole-to-bulge mass ratio K is fixed by the density 
at which the gaseous material in the inner region 
becomes self-gravitating and forms stars, and by 
the value of p CT it- However, unless loop orbits of 
sufficient eccentricity are created in the whole har- 
monic halo core region, the mass of the bulge will 



continue to grow. At the same time, the dynami- 
cal infall onto the SBH can choke, if the remaining 
gas mass is insufficient to cause bar instability and 
channel it to the SBH. In the latter case K will de- 
cline below its initial value as the bulge grows, but 
will tend to approach a final value roughly inde- 
pendent of Pcrit- In the opposite limit, the SBH 
will grow in proportion to the bulge, and K will 
stay constant (as shown by Figs. 2 and 4). In ei- 
ther limit, we assume that the process stops when 
the mass of the baryonic component reaches the 
minimal mass required to achieve loop orbits with 
a given p > p C rit, at all radii. 

For less triaxial haloes, and for a given bulge 
core density, lower SBH masses are required to 
produce loop orbits below a critical eccentricity. 
The minimal mass of the bulge needed to create 
the required loop orbits in the outer region, how- 
ever, also decreases — especially since, for a mildly 
triaxial halo, the radius at which no loop orbits 
exist inside the harmonic core decreases with de- 
creasing triaxiality. Thus, in the regime of mild 
halo triaxiality, the ratio of K to (normalized) 
minimal bulge mass should not depend sensitively 
on the value of (3. We do ignore the fact that tri- 
axiality can be a function of radius, and choose a 
fixed value of (3 for simplicity. 

3.5. Constraints on morphology and halo 
properties 

Once orbits of sufficiently large p are present, 
dissipation will reduce radial motion relative to 
these orbits, as well as vertical motion away from 
the symmetry plane defined by the angular mo- 
mentum vector. Infalling gas will start to popu- 
late the newly-formed, round, non-intersecting pe- 
riodic loops, leading to the formation of a disk 
component inside the halo core. At this stage 
the bulge formation stops. Outside the halo core, 
gas can accumulate at any stage of the formation 
process on closed loop orbits, which always exist 
in strongly inhomogeneous density distributions. 
Provided that the halo triaxiality is small, these 
orbits will be nearly circular (for potential axis 
ratio (3 mildly deviating from unity p ~ (3: see, 
e.g., Rix 1995). 

If the halo core radius is exceedingly large, how- 
ever, it is possible that no significant disk compo- 
nent will form at all. For if the core is large relative 
to the total halo (virial) radius, the contracting 
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gaseous component will end up in the core, instead 
of spinning up and forming an extended disk. Los- 
ing most of its angular momentum to the halo, it 
will eventually end up as part of the bulge compo- 
nent. This would lead to the formation of an ellip- 
tical rather than a disk-dominated galaxy. In this 
context, one expects that haloes with larger core 
radii are host to larger spheroidal components. 
Since it appears that, in general, more massive 
galaxies are usually of earlier type (e.g., Persic, 
Salucci & Stel 1996), one could deduce that more 
massive cores have larger core radii. This would be 
expected if these cores followed a Faber- Jackson 
type relation (for example, if Mh ~ cx^, then 
Mn ~ as tentatively suggested by observa- 
tions (Burkert 1995; Salucci & Burkert 2000; Dal- 
canton & Hogan 2001) and deduced if halo cores 
formed via the destruction of the inner (p <~ r _1 ) 
regions of NFW haloes (Navarro, Frenk & White 
1997). In this case, lower mass haloes form when 
the Universe is denser and are, therefore, more 
concentrated. As a consequence, the region where 
p r*j r _1 is smaller relative to the virial radius for 
low mass haloes (El-Zant, Shlosman & Hoffman 
2001). The existence of a Faber- Jackson type re- 
lation for halo cores is also required in order to 
reproduce the M. — ob relation as discussed in 
the next section. 

Within the above framework, the formation of 
the SBH is intimately tied to the bulge compo- 
nent, whereas the formation of the disk component 
takes place after the processes leading to bulge 
and SBH formation are essentially complete. The 
bulk of the disk is also expected to form at scales 
larger than the halo harmonic core (but see Sec- 
tion 5). Thus, the deduced correlation involving 
the bulge and SBH does not simply generalize to 
one involving the disk as well, in accordance with 
observations (e.g., Kormendy & Gebhardt 2001). 

4. The M. - <7 B relation 

The fact that all the orbital properties discussed 
in this paper depend only on the relative mag- 
nitudes of the masses and spatial scales of the 
galactic components involved has important con- 
sequences. For the M, — Mb relationship it has 
the obvious implication that the derived correla- 
tion will hold for all halo masses. A more powerful 
prediction transpires in relation to the M, — ctb 



relation. For this to hold, it is necessary that the 
masses and velocity dispersions of halo cores are 
related in a similar manner. 

Suppose that for a given halo there exist unique 
values for the SBH and bulge masses, and for the 
bulge scalelength. If a Faber- Jackson type relation 
between the halo core mass and velocity dispersion 
exists, a corresponding relation will exist between 
the bulge parameters. It also follows that a simi- 
lar relation will exist between the SBH mass and 
the bulge (and halo) velocity dispersion. This can 
be deduced by simple scaling transformations, be- 
cause the orbital properties we are interested in 
are all invariant with respect to spatial scale and 
mass transformations. This means that if we mul- 
tiply the masses of the SBH, bulge and halo by 
some constant o>m, thus effectively changing the 
mass units, the curves in Fig. 1 will remain invari- 
ant. In the same manner, if we multiply all length- 
scales by some factor «r, so that Ru ckr-Rh 
and i?B — * (xrRb, thus effectively changing the 
length unit, all curves in Fig. 1 remain the same. 
This implies that if, for example, Mh oc R^, then 
M. oc Mb oc Rq oc R^ are equivalent systems in 
the sense decribed above (they have the same axis 
ratio curves for their loop orbits with the same p 
values as a function of rescaled radius) . These sys- 
tems will all follow the Mh — oh relation for the 
halo. In this particular case, cth oc M 1 / 4 implies 
a similar relationship between ctb and Mb, as well 
as M.. 

It is of course possible that, for a given halo, 
the masses of the SBH and bulge, as well as the 
bulge lengthscale, are not unique. In other words, 
the subset of haloes with a given Mh and Ru may 
contain bulges and SBHs with a distribution of 
properties. In Section 3.2 we had assumed that 
the bulge collapse is largely terminated when its 
density is of the order of the halo core density. 
This in turn fixes K, once p cr it is determined. In 
this section we further assume that the value of K 
is not changed as a result of any subsequent bulge 
growth. That is, the SBH grows in tandem with 
the bulge (cf. Section 3.4). This relaxes the as- 
sumption that the bulge masses are determined by 
the minimal mass required to create closed loops 
with p > p cr it at all radii, allowing for more mas- 
sive bulges. 

One would like to infer to what extent varia- 
tions in the bulge and SBH properies, within a 
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given halo, affect the homology relations discussed 
above. We now show that when one accounts for 
variations of the bulge and the SBH parameters, 
the departure from the aforementioned relation is 
not dramatic. This comes about basically because 
the mass and core radius of the bulge are corre- 
lated, under the assumptions of our model. First, 
we note that the average velocity dispersion of 
the core of the baryonic component can be writ- 
ten as u\ = aGMB/i?B, and the average density 
PB = 2- 3 / 2 M B /i?|. In this case a B oc p B /6 M^ /3 . 
Here a depends on the functional form of the den- 
sity distribution. Its exact value is unimportant if 
all bulges are assumed to have the same functional 
form for their density distributions. Henceforth we 
set a = 1. 

For a constant bulgc-to-halo density ratio pb a 
relationship between bulge velocity dispersion and 
mass Mb oc <t b results. In a toy model where the 
halo density does not vary at all with radius a con- 
stant bulge density determines a unique K, given 
p cr it, and an "M — cr" relation between bulge prop- 
erties within a given halo arises, with index equal 
to 3. Considering more realistic models for the 
density distribution raises the index somewhat. 
When the core radius of the baryonic component 
is not very small compared to that of the halo, the 
density of the halo will not be strictly constant in 
the region of interest. As a result, the slope of the 
M — a relationship will further increase. 3 This oc- 
curs because larger bulges "see" a smaller mean 
halo density, and therefore smaller bulge densities 
suffice to produce the same relative contribution 
to the potential. Since Mb oc Ogpg 1 ^ ^ an m " 
verse correlation between density and a steepens 
the Mb — cb relation. (A marginal effect is already 
seen in Fig. 1 where the inner maxima produced 
by bulges with progressively larger Rb also have 
progressively larger inner minima of p.) 

If the baryonic component giving rise to the 
gravitational potential is significantly flattened 
there will be an increase in the M, — <tb slope 
to 4. In this case, the surface, rather than the vol- 
ume density, will determine the gravitational field, 
resulting in a situation where Mb oc <j b . 
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3 If pe °c R b ° , where a > 0, the slope of M. — o"b is given by 
2[(3 — «)/(2 — a)] which stays between 3 — 5 for a = — 1.5, 
and then rises rapidly. The resulting slope of M. — <tb is 
some weighted average of pb, and greater than 3. 



Fig. 7. — Relationship between bulge mass, Mb, 
and velocity dispersion, cr B , for systems with ps = 
1 and for different p C rit, as in Fig. 2. Open circles 
mark the lower ends of curves corresponding to dif- 
ferent values of p cr it, ranging from 0.3 to 0.8, left 
to right, and represent the increase of the minimal 
bulge mass with p crlt . The curves for different p cr it 
virtually coincide and the circles are placed to de- 
note the minimal mass, which is a distinguishing 
feature since it differs with p cr it- The dashed line 
has a slope of 4. The slope of the solid lines is 
about 3.2, which slightly increases for higher ve- 
locity dispersions. 

The fact that the index of the relation between 
Mb and <tb at constant bulge density is close 
to that of the Faber- Jackson relationship implies 
that, even if there is considerable variation in the 
bulge mass for a given halo, the departure from a 
Faber- Jackson relation defined by the halo param- 
eters (as described above) is not too large. This 
is illustrated in Fig. 7, where we plot the rela- 
tionship between the average velocity dispersion 
(simply defined as cr| = GMb/Rb) of the bary- 
onic component and its mass for systems having 
parameters K and p cri t corresponding to those in 
Fig. 2. The plots are obtained by keeping K and 
Pcrit constant and, for a given bulge mass, decreas- 
ing its characteristic radius until p > p cr ;t at all 
radii inside the halo core. To obtain physical pa- 
rameters we set Vh = 200 km s" 1 and Rb = 5 
kpc. 

The solid line in Fig. 7 is actually a superpo- 
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sition of several lines, corresponding to different 
values of p cr a (and K), as used in Fig. 2. Its aver- 
age slope is around 3.2 and increases slowly toward 
higher dispersion velocities. These lines, however, 
have different end-points, determined by the mini- 
mal bulge masses associated with the different p CI i t 
(cf. Section 3.3) which are denoted by circles in 
Fig 7. 4 The maximum values of the bulge masses 
in these curves correspond to 2.5 times the halo 
core mass. It seems unlikely that bulges would be 
more massive than this. Their characteristic den- 
sities or radii would have to be much larger than 
those of the halo, as would their contribution to 
the potential. In systems with disks such extreme 
conditions would violate the shapes of observed 
rotation curves. 

In reality, for a given halo, there is a range of 
possible bulge densities and associated values of 
the K parameter. This, in principle, affects the 
normalization of the Mb — ob relation. Never- 
theless, the normalization turns out to be weakly 
dependent on p B (~ Ob/Ob^ 2 ). Thus, variations 
in bulge densities within haloes of given mass and 
size do not significantly affect the relationship be- 
tween bulge mass and velocity dispersion. In our 
model, the central (bulge) density required for the 
production of loop orbits with p > p„it in the scal- 
ing regime is also only weakly dependent on K (as 
K 2 / 5 , cf. Fig. 2). The above leads to an important 
corollary, that large variation in K will cause only 
small changes in Mb — <tb, as illustrated in Fig. 8. 
This relation, therefore, appears to be robust and 
is not heavily affected by changes in bulge and 
SBH parameters. Note also that, for less dense 
bulges the slope of the lines in Fig. 8 tends to ~ 4, 
especially in the limit of large values of the velocity 
dispersion (which is also where most of the obser- 
vations lie). This is due to the the effect described 
above — massive bulges with smaller densities are 
more extended, and, therefore probe larger regions 
of the halo core. 

To obtain an M # — ob relation, one has to multi- 
ply the values of the bulge masses in Fig. 7 by the 

4 We note here that when the minimal mass is approached, 
it is necessary to contract the bulge by a large amount to 
maintain p > p cr it as the mass decreases. The result is 
that small "hooks" (not shown here) appear in place of the 
circles in Fig. 7. In principle, this regime may be of interest, 
but is not considered in this paper because it occurs for a 
very small range of bulge masses. 
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Fig. 8. — Relationship between bulge mass, Mb, 
and velocity dispersion, ctb, for systems with 
Peru = 0.8 but with different K (and, therefore, 
different densities, as would be clear by following 
the line with p cr j t = 0.8 in the scaling regime of 
Fig. 2). Solid lines, from bottom to the top, are 
K = 1CP 6 increasing to 1CP 2 . The left /right ex- 
tensions of these lines are defined by the limits 
on Mb (see text). The dashed line has a slope of 
4. The slope of the solid lines is variable, ranging 
from ~ 3 for small a and K, to ~ 4 for large values 
of these parameters. Note the small dispersion in 
Ob for large variations in Mb and K. 

appropriate K factors, arriving at an M. — ctb re- 
lation within a given halo. This introduces "scat- 
ter" by shifting the K — const, lines, as can be 
seen from Fig. 9. For any fixed value of ctb, this 
figure reveals that the (vertical) scatter in M. is 
similar to that in K . 

This leads us to an important point. In the 
present formulation, the Mb — ob relation is 
tighter than the M. — ob one. This appears to 
contradict observations suggesting that M. — er B 
is much tighter than the Faber- Jackson relation. 
However, this result has been obtained under the 
assumption that there is no scatter in the halo 
Faber-. Jackson relation. A significant scatter in 
Mh — oh would result in a corresponding scat- 
ter in M B — ctb ! which follows from the homology 
scaling discussed in the beginning of this section. 
Moreover, if p cr i t correlates with cr B in such a way 
that gaseous sustems embedded in haloes with 
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Fig. 9. — Same as in Fig. 7 but for a range of 
M,, assuming zero scatter in the halo's Mr — oh 
relation. The solid lines are, therefore, shifted ac- 
cording to their K values, increasing from bottom 
to the top. These correspond to K = 5.7 x 10~ 5 , 
1.2 x 1CT 4 , 2.6 x 10 4 , 5.9 x 1(T 4 , 1.6 x 1CT 3 and 
K = 6.2 x 10~ 3 (as can be inferred from Fig 7, 
by considering the intersections of the lines with 
different p crit with the dashed line representing 

PB = 1). 

higher velocity dispersion require larger values of 
p„it to remain stable, the scatter in the M, — ctb 
relation due to variations in K can be reduced. For 
example, in Fig. 9, we have assumed a single halo 
for all the lines with different p a it- If, however, 
lines with higher values of p cri t are associated with 
haloes with larger cth (and thus ctb), the constant 
p cr it lines would be shifted in such a way that the 
vertical distances between them at a given value 
of Ob decrease. We illustrate this effect in Fig 10, 
where a linear relationship is assumed between ctb 
and p crit - 

We suggest that the relative tightness of M, — 
ctb could result from a loose relation between ct b 
and p cr it, coupled with scatter in M H — ct h com- 
parable to that of the Faber-Jackson relation for 
bulges. Such p cr it — ctb relation can result from a 
general trend that dissipation is increasing in more 
concentrated systems. Testing this will require de- 
tailed modeling of the way in which the dissipation 
rate depends on p cr n, velocity and density, and is 
outside the scope of this paper. Based on the il- 



Fig. 10. — Reducing scatter in M. — ctb by 
Pcrit — ctb correlation. A linear correlation between 
Pcrit and &b was imposed in Fig. 9 (in practice 
this is done by multiplying each value of ctb by 
Pcrit/0.5). This in effect assumes scatter in the 
halo Faber-Jackson relation, implying that, for 
any given mass, a range in halo (and, therefore, 
bulge) velocity dispersion is plausible and that sys- 
tems embedded in more concentrated haloes are 
associated with higher values of p C rit- The solid 
lines are, therefore, shifted according to their val- 
ues of p cr it , significantly reducing the vertical dis- 
tance between them. 

lustrative example of Fig. 10, we only claim here 
a plausibility of a loose relationship between p cr j t 
and the system velocity dispersion. If the pos- 
tulated correlation persists for haloes of different 
masses, a steepening of the M. — ctb relation rela- 
tive to the halo (and bulge) mass velocity disper- 
sion relation is expected, implying that the former 
should have a larger index. 

5. Discussion and conclusions 

In the model presented here, gaseous baryonic 
material settles inside a mildly non-axisymmctric 
halo with a nearly constant density core. Initially, 
no orbits with a definite sense of rotation exist. 5 
The first infalling baryonic material, therefore, ef- 
ficiently loses its angular momentum to the core. 
This initial collapse terminates only when the gas 

'Or only very eccentric ones if one allows for slow figure 
rotation in the halo. 



14 



becomes self-gravitating and forms stars. The in- 
creased central density concentration which is pro- 
duced in this first phase, however, destroys the 
harmonic core, paving the way for the existence of 
non-intersecting closed loops with a definite sense 
of rotation. 

If the loop orbits are too eccentric, the gas will 
shock and depopulate them. These orbits, there- 
fore, cannot represent long-lived attractors of the 
dissipativc motion. Thus, if the baryonic com- 
ponent does not possess a central cusp initially, 
the absence of sufficiently round supporting closed 
loops will lead to the formation of a central mass 
concentration including an SBH. At larger radii 
the potential responsible for creating sufficiently 
round loop orbits is that of the extended baryonic 
component, in the form of a bulge. This leads to 
a linear relationship between the bulge and SBH 
masses. 

It is crucial that both the onset of self-gravity 
and the appearance of increasingly circular loop 
orbits are subject to the same condition, both de- 
pending on the density ratio of the collapsing gas 
to that of the background halo core. This lim- 
its the allowable range in densities. For bulge 
core densities of order those of the halo core, and 
for plausible values of the critical eccentricity, the 
value of K lies in the range 1CP 4 — 10~ 2 , which is 
compatible with present observations. For densi- 
ties up to an order of magnitude larger, the range 
of values is the same, provided that the critical val- 
ues of the closed loop axis ratios are larger. This 
can be expected if the dissipation rate along loop 
orbits is dependent on density — a plausible as- 
sumption. 

The most important prediction of the model 
outlined above is that the bulge and SBH param- 
eters are determined by the halo properties. In 
particular, relationships between SBH and bulge 
masses and the velocity dispersion of the bulge 
necessarily arise, and with the right exponents, 
if the haloes should also exhibit a Faber-Jackson 
type relation between their masses and velocity 
dispersions. Moreover, within a given halo, vari- 
ations in the bulge and SBH properties are not 
expected to destroy these relations — because 
imposing the condition of critical eccentricity re- 
quires that the bulge and SBH masses are related 
to the bulge velocity dispersion via a power law, 
with index also close to that of the Faber-Jackson 



relationship. Finally, if the Faber-Jackson rela- 
tionship for the halo exhibits significant scatter 
and if, as again seems plausible, the critical eccen- 
tricities of the loops anti-correlate with the den- 
sity of the system, less scatter should be present in 
the relationships between SBH masses and bulge 
velocity dispersions than the corresponding rela- 
tionships between bulges masses and their own ve- 
locity dispersions — the standard Faber-Jackson 
relationship. 

There is already tentative evidence that halo 
cores may indeed follow Faber-Jackson type re- 
lationships (Burkcrt 1995; Dalcanton & Hogan 
2001). In addition, a halo core produced by flat- 
tening out the inner region (where p oc r" 1 ) of 
the NFW profile would produce such a relation 
(El-Zant, Shlosman & Hoffman 2001). The model 
thus makes testable predictions concerning the re- 
lationships among the very inner regions of galax- 
ies, their extended baryonic components, the dark 
matter haloes they are thought to be embedded in 
and the cosmology that predicts their existence. 

In this framework, larger and more massive halo 
cores produce, on average, larger and more mas- 
sive bulges. Disks form outside the core, or later, 
when the central concentration produced by the 
baryons destroys the core. If the core is very large, 
most of the baryonic material is consumed in the 
first phase and no significant disk forms. This ef- 
fect is expected to be prominent in larger mass 
cores, since if these follow the Faber-Jackson rela- 
tion, more massive haloes should have proportion- 
ally larger cores. Other predictions include the 
requirement that the average density of the bulge 
in the central region should be close to that of the 
halo core. There is also a minimal bulge mass as- 
sociated with a given core, although this varies sig- 
nificantly with the critical loop orbit eccentricity 
assumed. A number of additional consequences for 
galaxy formation and evolution will be discussed 
elsewhere. 

To obtain relationships between SBH masses 
and bulge properties, we have assumed that dur- 
ing the gaseous infall phase the baryonic compo- 
nent did not have a central density cusp. Most 
observed bulges and ellipticals, however, do have 
such cusps. One then has to assume that once 
star formation starts, cold dissipational collapse 
is initiated. This would lead to a central density 
cusp, as it does in the case of cosmological haloes. 
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Memory of the initial state is retained via the total 
energy, which determines the final velocity disper- 
sion. 

Merritt & Ferrarese (2001) have suggested 
that a "self-regulating" mechanism, related to the 
threshold mass necessary for the loss of triaxiality 
in the system, may be behind the close correlation 
between SBH and bulge properties. As these au- 
thors point out, however, the SBH masses required 
for strong chaotic behavior leading to rapid loss of 
triaxiality are in fact probably too large — being 
of the order of a few percent of the mass of the 
system's baryonic component. This is actually of 
the order of the SBH mass needed to create round 
loop orbits at all radii inside the halo core. In 
our model, however, this is not assumed. Instead, 
an additional baryonic component plays the role 
of creating these orbits in the outer region. The 
collusion between this "bulge" component and the 
SBH in destroying the harmonic core and creating 
a situation whereby stable gaseous motion can ex- 
ist gives rise to the correlations described in this 
paper. 

The correlations obtained here are compatible 
with the observed ones, and with acceptable scat- 
ter, despite our lack of knowledge of the values 
of such parameters as p cr it and its variation with 
system properties. Detailed modeling of the gas 
dynamics will be required to further constrain this 
model. It is also possible that our distinction be- 
tween the dynamical role played by the SBH and 
that played by the bulge core is too restrictive. We 
have introduced this to be able to obtain quanti- 
tative results, within the model, solely on the ba- 
sis of the orbital characterisitics. In general, the 
roles of the two components, i.e., the SBH and the 
bulge, may not be too distinct — formation of the 
SBH can take place simultaneously with a cuspy 
bulge. For this, the even more ambitious task 
of a self-consistent treatment, including gas and 
stellar dynamics and star formation, is required. 
We believe, however, that our results are generic 
and arise from fundamental dynamical phenomena 
which will manifest themselves in any formulation 
of galaxy formation in mildly triaxial haloes with 
harmonic cores. 
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A. Perturbation analysis 

Let the potential $ be a function of the Cartesian coordinates in the halo symmetry plane, $ = $>(x 2 , y 2 ) 
with $(0,0) = 0. We follow de Zeeuw & Merritt (1983, hereafter dZM) and expand the potential in even 
powers of the coordinates: 

* = \k\x 2 + l;K 2 2 y 2 + ja 5 x 4 + ^a 7 x 2 y 2 + ja g y 4 + ... (Al) 

If no bulge exists, inside the halo's harmonic core only the terms quadratic in the coordinates are important. 
If, due to an additional component, or at the boundaries of the core, weak nonlinearity is present the series 
can be truncated, as above, to second order, as higher order terms are unimportant. The coefficients arc 
given by: n\ = 2d x 2 k 2 — 2<9 y 2$, a$ — 2<9 2 2 $, 07 = 2d 2 . 2y 2$, ag = 2d 2 2 ^, where the derivatives are taken 
at x = y = 0. In terms of these one defines the auxilliary variables /in = §^|, M12 = 5^7^; ^22 = f fir- 
The condition for stability of the loop orbits is (sec dZM; Table 2A third row) 

Mi2(Mn - M12 + M22) > 0. (A2) 

We will be interested in systems that are both mildly nonlinear and mildly nonaxisymmetric: thus K\ <~ n 2 
and as ~ ag <~ aj < 0. The above condition reduces to 

3 f n 2 Ki \ , . „. 

2 \m n 2 J 

which, under the above conditions, is always satisfied. 

It is, therefore, the condition for the existence of loop orbits and not their stability that will be of interest 
to us. In the absence of a central mass (SBH) or density cusp these cannot be found arbitrarily close to 
the center, instead there is a bifurcation radius beyond which these exist. This determines the effective core 
radius of the system. To second order, the condition for the existence of loop orbits is 

1 k 2 S 1 

M22 - 2^12 -~q - 2^ 12 ~ ( ) 

with 5 — ki/k 2 — 1 and (to first order) Q = H/n 2 , where H is the Hamiltonian. For H > (first order 
potential terms dominate in eq. Al) and n\ < k 2 , the above requires 

H > 4 "i"lte-"i). (A5) 
n 2 a7 — 3niag 

In the unperturbed case, the action variables (e.g., Binney & Tremaine 1987) are time-independent and 
(exact) solutions can be written in terms of these in Cartesian coordinates as x — yj2I\/ni cos#i, etc. In 
the mildly nonlinear case the solutions of the equations of motion, averaged over a dynamical time (denoted 
below by a "bar"), approximate the true solution to first order in the (relative amplitude of the) perturbation 
and for a number of dynamical times inversely proportional to this (dZM; see also Bogoliubov & Mitropolsky 
1961; Arnold 1989). In this case analogous approximate solutions can be given in terms of the corresponding 
action variables. In particular, for loop orbits in mildly nonlinear potentials one finds 

f Q(|Ml2 - M22) + Kl - K 2 

h = ; (Ab) 

-Mil + M12 - M22 

and 

h = Q-h- (A7) 

Eliminating Q one gets 
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where f\ = -/in + fii 2 - M22, F 2 = «i — k 2 and F 3 = fi i2 - fi 22 . 
At bifurcation, to first order, 

h = Q = H/n 2 , (A9) 

where is given by eq. (A5). These orbits are infinitely thin and represent oscillations along the y-axis, 
with amplitude 

2W = ^2I 2 /k 2 , (A10) 

which is the effective core of the bulge-halo system, as mentioned above. As one increases H these become 

thicker with axis ratio x max /y max = The view we have taken in this paper is that when this ratio 

becomes large enough, > p CT u, such loop orbits can support gaseous motion, and that stars formed on these 
orbits can constitute populations of stellar disks, thus ending the bulge formation stage. We note here that, 
beyond the bifurcation point, the above relation predicts a rather rapid (as a function of radius) transition 
to round loop orbits. This is confirmed by orbital integration, even though if the bulge mass is smaller than 
a certain minimal mass, in the fully nonlinear treatement, the orbital axis ratio can decrease again. The 
second minimum in axis ratio curves (e.g., Fig. 1) is thus not reproduced by this perturbation analysis. 



B. Application to the potential used in this paper 

We now apply the perturbation analysis to the superposition of potentials given by eqs. (1) and (3) with 
the goal of calculating the bifurcation radii of the loops. This is the effective harmonic core radius of the 
bulge-halo system. Since we are only interested in orbit shapes, we set V B = G = 1 and use scaled variables, 
as defined in Section 2. For the halo, we use 

$ H = ll g(l + x 2 + qy 2 ), (Bl) 

where q, corresponding to /3~ 2 in eq. (1), parametrizes the nonaxisymmetry. For the bulge, we have 

*b = - , h ■ (B2) 



With the above definitions for the potential we have 



K 2 = 2d x ^ = l + ^-, (B3) 
4 = 2d y ^ = q+^-, (B4) 



a 7 = 2d 2 x2y ^ = -q- 3 -^, (Br,, 

and 



a 9 = 2d 2 2 ^ = - q i--^L. (B6) 

From equations (A5), (A9), and (A10) bifurcation happens when the value of the long axis of the loops 
satisfies 

2 = 8kik 2 (k 2 - Kl) 
^ max K 2 a 7 -3K ia9 ' 1 ' 
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Substituting from the expressions above and taking the limits Mb/R%{~ Pb) 3> 1, Mb/R% ^ 1> an d 
assuming q ~ 1, we obtain 

,, 2( g -l)V» i% 2 
W-^g— (B8) 

This is the effective harmonic core radius of the bulge-halo system. 
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